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Abstract 



In this talk I report on the results of a recent calculation of the corrections 
to a three-point correlation function at the three-loop level in QCD perturbation 
, theory, which allows one to extract the matrix element of — mixing with 

next-to-leading order accuracy pQ. The evaluation of mixing parameter at NLO 
fS| ' allows for a consistent analysis of B^ — B^ mixing since the coefficient functions of 

. the effective Hamiltonian for this process are known with the necessary accuracy. 

Q ■ Presently the pattern of CP violation in the Standard Model is under thorough experi- 

^ . mental study in dedicated experiments by BABAR collaboration at SLAC (e.g. PlEl) 

Q-i! and BELLE collaboration at KEK (e.g. [HE])- The standard mechanism of CP violation 

{iLii to be primarily tested is a Cabibbo-Kobayashi-Maskawa paradigm with mixing of (at 

least) three quark generations (for review, see e.g. P). At the hadronic level the fact of 
quark mixing mainly reveals itself as mixing of neutral pseudoscalar mesons. The most 
. famous system where mixing occurs and has been studied in much detail is the system 

\ of neutral kaons. Study of K'^ — mixing strongly constrained the physics of heavy 

\ particles and allowed to estimate the numerical value of the charm quark mass from 

the requirement of GIM cancellation before the experimental discovery of charm (see 
e-g- [ZD- At present the experimental studies of CP violation shifted to the realm of 
heavy mesons for which they are considered more promising. In particular, recent ex- 
perimental results for heavy charmed mesons D{uc) are encouraging [8j. However the 
systems of B(i{dh) and Bs{sb) mesons are the most promising laboratory for performing a 
precision analysis of CP violation and mixing both experimentally and theoretically jU] . 

Mixing in any system of neutral pseudoscalar mesons is described by a 2x2 effective 
Hamiltonian or mass operator (M — ir/2)^j, {i,j} = {1,2} where M is related to the 
mass spectrum of the system and F describes the widths of the mesons. In the presence 
of flavor violating interactions {AB = 2 in our particular case) the effective Hamiltonian 
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acquires non-diagonal terms. The difference between the values of the mass eigenstates 
of B mesons Am = M^eavy — Mug^t ~ 2 \Mi2\ is precisely measured Am = 0.489 ± 
0.005{stat) ± 0.007{syst) ps~^ jTU]. With an accurate theoretical description of the 
mixing, it can be used to extract the top quark CKM parameters. In Wolfenstein's 
parametrization the CKM matrix elements reveal a hierarchy in magnitude. In terms 
of Cabibbo angle A ~ 0.22, this parametrization of the CKM matrix V reads 



V 



( 1 - |A2 
\A\^ {\ - p - iri) 




Here 



A, F,fe = AA^ V^^ = A\\p-%r]), Vu = A\\\-p-if]). 



While |V^b| and |\4b| can be extracted from semileptonic B decays, \ytd\ is at present 
probed in the process of B\-B\ mixing. 

The expression for the effective Hamiltonian describing AB = 2 transitions is known 
at next-to-leading order (NLO) in QCD perturbation theory of the Standard Model [TT] 
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where Gp is a Fermi constant, Mw is the W-boson mass, rjB = 0.55 ± 0.1 [12], J5 = 
1.627 in the naive dimensional regularization (NDR) scheme, S'o(x() is the Inami-Lim 
function [12], and = (&L7(T'^L)(&L7o-c?L)(Ai) is a local four-quark operator at the 

normalization point p. Mass splitting of heavy and light mass eigenstates is 



Am = 2\{B'\H^^='\B')\=C[af\p) 
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where we introduced a constant C = GpM^, {Vtb*Vtd)^ rjBmBSoixt) / (47r^). The largest 
uncertainty in calculation of the mass splitting is introduced by the hadronic matrix 
element A = (-B°|0(/x)|5°) that is poorly known JU]. The evaluation of this matrix 
element is a genuine non-perturbative task, which should be approached with some 
non-direct techniques. The simplest approach ("factorization") ^3] reduces the matrix 
element A to the product of simpler matrix elements measured in leptonic B decays 



M = ^{B%bLl.dL\0){0\bL7''dL\B') = -flm\ 



(1) 



where the decay constant fs is defined by (0|6L7^(ii|i?°(p)) = ip^fB/2 and is the 
B^ meson mass. A deviation from the factorization ansatz is usually described by the 
parameter Bb defined as ^ = BbA-^; in factorization Bb = ^- The evaluation of this 
parameter (and the analogous parameter Bk of K^ — K^ mixing) has long history. Many 
different results were obtained within approaches based on quark models, unitarity. 
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ChPT. The approach of direct numerical evaluation on the lattice has also been used. 
The corresponding results can be found in the literature [1^1 HSl HZl UHl UH 1201 1211 122] • 

In my talk I report on the results of the calculation of the hadronic mixing matrix 
elements using Operator Product Expansion (OPE) and QCD sum rule techniques for 
three-point functions [H UHl HZl (HI I2H1 121] • This approach is very close in spirit to 
lattice computations j2I], which is a model- independent, first-principles method. The 
difference is that the QCD sum rule approach uses an asymptotic expansions of a Green's 
function computed analytically while on the lattice the function itself can be numerically 
computed provided the accuracy of the technique is sufficient. The sum rule techniques 
also provide a consistent way of taking into account perturbative corrections to matrix 
elements which is needed to restore the RG invariance of physical observables usually 
violated in the factorization approximation 

To start with let me introduce the three-point correlation function 

n(j9i,p2) = / dxdy{0\TJs{x)O{0)JB{yme'''''-''^' (2) 

of the relevant AB = 2 operator 0{fi) and interpolating currents for the i?°-meson 
Jb = (jrib + mdjdi'y^b. Here nib is the b quark mass. The current Jb is RG invariant and 
Jb = d^i^d-y^-y^b). The main relevant property of this current is {0\Jb{0)\B^{p)) = fBfn% 
where is the 5- meson mass. A dispersive representation of the correlator reads 

\—-nf 2 2 2\ f p{.Si,S2,q^)dsids2 , . 

Uip,,p2) = Uip„p2,q ) = J (^^_^.)(^^_^.) (3) 

where q = P2 — Pi- For the analysis of — B^ mixing within the sum rule framework 
this correlator can be computed at = 0. 

Phenomenologically the matrix element {B^\0{fi)\B^) determines the contribution 
of the 5-mesons in the form of a double pole to the three-point correlator 

Wi,pI = + . . . (4) 

m% -pi m%- pi 

Because of technical difficulties of calculation, a practical way of extracting the B^ — B^ 
matrix element is to analyze the moments of the correlation function at = = at 
the point = 

d'+m{pl,pl,0) f p{si,S2,0)dsids2 



A theoretical computation of these moments reduces to an evaluation of single scale vac- 
uum diagrams and can be done analytically with available tools for the automatic com- 
putation of multi-loop diagrams. Note that masses of light quarks are small (e.g. |26j ) 
and can be accounted for as small perturbation which is relevant for the problem of 
Bg meson mixing 127]. The leading contribution to the asymptotic expansion is given 
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Figure 1: Perturbation theory diagram at LO 

by the diagram shown in Fig. ^ At the leading order in QCD perturbation theory the 
three-point function Il{pi,p2) of Eq. © completely f actor izes 

n(pi,p2) = ^uMw^ip^) (5) 

into a product of the two-point correlators n^(p) 

U^{p) = p^U{p') = J dxe'P^{0\TJs{x)hi^dLmO). (6) 

At LO the calculation of moments is straightforward since the double spectral density 
p{si, S2,q'^) is explicitly known in this approximation. Indeed, using dispersion relation 
for the two-point correlator 

one obtains the LO double spectral density in a factorized form 

8 4 
P^"(si, S2, q^) = 3 (Pi ■ P2)p{si)p{s2) = -{si + S2- q^)p{si)p{s2) . (8) 

Thus, all PT contributions are of the factorizable form at LO. First non-factorizable 
contributions to Eq. Q appear at NLO. Of course, at NLO there are also the factorizable 
diagrams. Note that the classification of diagrams in terms of their factorizability is 
consistent as both classes are independently gauge and RG invariant. 

Consider first the NLO factorizable contributions that are given by the product of 
two-point correlation functions from Eq. as shown in Fig. |21 Analytical expression 
for such contributions can be obtained as follows. Writing n(p^) = IIloIp^) + nNLo(p^) 
one finds 

n^Lo(Pl.P2) = ^(Pl-P2)(nL0(p?)nNL0(P2) + IlNL0{pi)^L0{pl)). (9) 

Since the spectral density of the correlator HnloIp^) is known analytically the problem 
of the NLO analysis in factorization is completely solved. Even a NNLO analysis of 
factorizable diagrams is possible as several moments of two-point correlators are known 
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Figure 2: Factor izable diagrams at NLO 

analytically. Others can be obtained numerically from the approximate spectral den- 
sity |28j. 

The NLO analysis of non-factorizable contributions within perturbation theory is the 
main point of my talk. The analysis amounts to the calculation of a set of three-loop 
diagrams (a typical diagram is presented in Fig.|2)). These diagrams have been computed 
using the package MATAD for automatic calculation of Feynman diagrams [21] ■ The 
package is applicable only for computation of scalar integrals, the decomposition of 
the three-point amplitude into scalars is known jHOj- A scalar function ^4(^1,^2) of two 
four- momenta pi and p2 can be expanded in a series over the scalar variable pf, p\ and 
Pi ■ P2 of the general form 

A{puP2)= E (^n,n,Mriply'iPi-P2y'. 

jl,j2,j3 

For the coefficients ajijajs finds 



1 (-1)"-^^ 



2n2nj^\ 2^>'+^^-n(2j, + - 2k)\ 

minik,k+n~n) (2ji+j3-2A; + 20! 
, . . - A; + /)! 

l=max(0,k—ji) ^-^ ^ ' 

(1 - £)2ji+j:i-2fc+2/+l(l - ^)ji+j3~k+l 
(1 - £)2ji+j.s-2k+2l{l - 2e)2j^+j3-2k+2l+l 

(1 - £)2ji+j3-2fc+;(l - - 2e)i 

(1 - £)2ji+j3-fc+«+l(l - £)jl+j2+j3-k+l+l 

X^11^22 ^12 ^[PliP2)\pi=p2=0 



with 



2 



& 
(^PiP'i 



Here (a)^ stands for the Pochhammer symbol (a)i, = r(a + 6)/r(a) and c? = 4 — 2e is 
the spacetime dimension. The actual calculation has been performed with the computer 
algebra system FORM I3T]. 
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Figure 3: A typical non-factorizable diagram at NLO 

Renormalization of the local four-quark operator O entering the effective Hamil- 
tonian has been done in dimensional regularization with an anticommuting 75. The 
renormalization of the operator O reads 

= _ ^1(9, (10) 

with C»c = {bLTf,^at''dL){bLT^"'''t''dL). The color space matrices are the SU^iS) gen- 
erators and T^^a = {ifilula — lalvln)/'^- Note that in four dimensional space-time the 
matrix T^^a reduces to the expression —ie^uapl^l^- This relation is however ill defined 
in Z)-dimensional spacetime of dimensional regularization. The renormalization of the 
factorizable contributions reduces to that of the 6-quark mass m. We use the quark pole 
mass as a mass parameter of the calculation. 

The expression for the "theoretical" moments reads 

where the quantities ajj, blj and Wj represent LO, NLO factorizable and NLO nonfac- 
torizable contributions as shown in Figs. ^El The NLO nonfactorizable contributions 
b^- with i + j < 7 are analytically calculated in ref. [Ij for the first time. The calculation 
required about 24 hours of computing time on a dual-CPU 2 GHz Intel Xeon machine. 
The calculation of higher moments is feasible but requires considerable optimization of 
the code. This work is in progress. The analytical result for the lowest finite moment 
Mt;,(2,2) reads 




^ 8366187 , 84608 233197 426319 

O99 = O2 G TT . 

17500 ^ 875 52500 315000 

Here S2 = ^^Ch (f) = 0.2604. . ., (3 = C(3), and //^ = m^. Higher moments contain 
the same transcendental entries S2, (3, vr^ with different numerical coefficients. The 
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Figure 4: A model of phenomenological spectrum 



numerical values for the moments are h"^- : ^2(^345) = {0.68, 1.22, 1.44, 1.56} and ^3(^4) = 
{1.96,2.25}. The above theoretical results are used to extract the non-perturbative 
parameter Bb from the sum rules analysis. 

The "phenomenological" side of the sum rules is given by the moments which can 
be inferred from Eq. (j3]) 

M ,)-^-FI I .0 ( 1 1 \ 

IVW,J) 3^i^^g.+,) +^[^ml + AY^Wi ' (ml + Ay^^mj^ ) 

where the contribution of the 5-meson is displayed explicitly. The remaining parts are 
the contributions due to higher resonances and the continuum which are suppressed due 
to the mass gap A in the spectrum model. A rough picture of the phenomenological 
spectrum is given in Fig. (jl)). 

For comparison we consider the factorizable approximation for both "theoretical" 

and "phenomenological" moments, which, by construction, are built from the moments 
of the two-point function of Eq. (0) 

According to standard QCD sum rule technique, the "theoretical" calculation is dual to 
the "phenomenological" one. Thus, Eq. (fT^ should be equivalent (in the sum rule sense) 
to Eq. (fTT|) . Also, in factorization, Eq. (fT^ is equivalent to Eq. (fT^ . Now Eq. (fT^ 
and Eq. (jlip differ only due to non-factorizable corrections. Therefore, the difference 
between Eq. (fT^ and Eq. ()13|) is because the residues differ from their factorized values. 
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To find the nonfactorizable addition to Bb from tlie sum rules we form ratios of the 
total and factorizable contributions. On the "theoretical" side one finds 

, '-^^ = 1 + — '2 — ^. (15) 

This ratio is mass-independent. On the "phenomenological" side we have 

Mphii, j) Bb + Rb{z' + z') + Cbz'+^ 



mIS, j) 1 + R^iz^ + z') + Cfz'+'J 



(16) 



where z = m% / {m?B + A) is a parameter that describes the suppression of higher state 
contributions. A is a gap between the squared masses of the 5-meson and higher states. 
Rb. Cb, R^ and are parameters of the model for higher state contributions within the 
sum rule approach. In order to extract the non-factorizable contribution to 5^ we write 
Bb = 1 + Ai?. Similarly, one can parameterize contributions to "phenomenological" 
moments due to higher 5-meson states by writing Rb = R^ + AR and Cb = + AC 
Clearly, AB = AR = AC = in factorization. We obtain 

Mph{i,j) _ Ai? + ARjz^ + z') + ACz'^^ 

mIS.J) l + Rf{zi + z') + CSz'+i ■ ^ ' 

Comparing Eqs. (|15|) and (jl7p one sees how the perturbative non-factorizable correction 
h^- is "distributed" among the phenomenological parameters of the spectrum. We ex- 
tract AB by a combined fit of several "theoretical" and "phenomenological" moments. 
The final formula for the determination of AB reads 



^hli = AB + AR(z'~^ + z'-^) + ACz'+^-^ (18) 



where AR and AC are free parameters of the fit. We take A = O.Am'^^ for the B meson 
two-point correlator. This corresponds to the duality interval of 1 GeV in energy scale 
for the analysis based on finite energy sum rules |32] ■ The actual value of AB has been 
extracted using the least-square fit of all available moments. Estimating all uncertainties 
we finally find the NLO non-factorizable QCD corrections to AB due to perturbative 
contributions to the sum rules to be 

A5 = (6±l)"^("^) 



An 

We checked the stability of the sum rules which lead to a prediction of AB. For = 
4.8 GeV and ^^(mfc) = 0.2 [33] one finds AB = 0.1. The calculation can be further 
improved with the evaluation of higher moments. The result is sensitive to the parameter 
z or to the magnitude of the mass gap A used in the parametrization of the spectrum. 

In conclusion, the B^ — B^ mixing matrix element has been evaluated in the frame- 
work of QCD sum rules for three-point functions at NLO in perturbative QCD. The effect 
of radiative corrections on 5^ is under complete control within pQCD and amounts to 
approximately +10% of the factorized value. 
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